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ABSTRACT 

We derive estimates for the cosmological bulk flow from the SFI++ Tully-Fisher 
(TF) catalog. For a sphere of radius 40/i~^ Mpc centered on the MW, we derive a bulk 
flow of 333±38 kms"^ towards Galactic {I, b) = (276°, 14°) within a 3° la error. Within 
a WOh-^ Mpc we get 257±44 kms"^ towards {I, b) = (279°, 10°) within a 6° error. These 
directions are at a 40° with the Supergalactic plane, close to the apex of the motion of 
the Local Group of galaxies after the Virgocentric infall correction. Our findings are con- 
sistent with the ACDM model with the latest WMAP best fit cosmological parameters. 
But the bulk fiow allows independent constraints. For WMAP inferred Hubble param- 
eter h = 0.71 and baryonic mean density parameter Qf, = 0.0449, the constraint from 
the bulk fiow on the matter density Qrm the normalization of the density fiuctuations, 
as, and the growth index, 7, can be expressed as as^Tn^'^^ = 0.86 ±0.11 
(for nm « 0.266). Fixing as = 0.8 and Qm = 0.266 as favored by WMAP, we get 
7 = 0.495 lb 0.096. The constraint derived here rules out popular DGP models at more 
than the 99% confidence level. Our results are based on a method termed ASCE (All 
Space Constrained Estimate) which reconstructs the bulk flow from an all space three 
dimensional peculiar velocity field constrained to match the TF measurements. At 
large distances ASCE generates a robust bulk flow from the SFI++ that is insensitive 
to the assumed prior. For comparison, a standard straightforward maximum likelihood 
estimate leads to very similar results. 

Subject headings: Cosmology: large-scale structure of the Universe, dark matter, cos- 
mological parameters 
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Introduction 



Cosmological bulk flows are the peculiar velocities of whole spherical regions around us. Bulk 
flows are usually considered for sufficiently large spheres where linear expressions for the velocity 
and density power spectra are valid. This greatly facilitates the calculation of expected bulk flows in 
cosmological mo dels, in contrast to analyzing the full field field which may involve non-linear effects 



on sm all scales (jFeldman et al.ll20ld : lAbate Erdogdull2009l : IZaroubi et al.ll200ll : iFreudhng et al 



19991 ) . In linear theory, the bulk flow of a sphere is solely determined by the gravitational pull of only 
the dipole component of the external mass distribution. Bulk flows are, therefore, an unmistakable 
indicator of distant large mass concentrations should they exist. The exact expression of the bulk 
flow, B(r), of a sphere of radius r is, 

3 



B(r) 



v{x)d X . 



(1) 



x<r 



where v{x) is the 3D peculiar velocity field as a function of the comoving coordinate x. Beneath this 
innocuous expression lie a multitude of nuisances. An unbiased estimate of B requires knowledge 
of V sampled uniformly overall the volume. However, observational probes of peculiar velocities 
measurements are available only for a few thousand galaxies with a patchy coverage of the local 
Universe. Further, peculiar velocity probes such as the TF relation allow us to constrain only the 
radial component of the peculiar velocities of galaxies. 

Recently compiled data on peculiar vel ocities has triggered renewed interest in the analysis 



of lar g e scale flows, includin g the bulk flow (jPavis et al.l l201ll : iLavaux et al.l l20ld : lErdogdu et al 
20061 ). iFeldman et al.l ( 2010l ) report an unusually large bulk flow of 416 ± 78 kms~^ in a sphere of 
100h~^ Mpc which is at odds with the ACDM model with the best fit pa,rameters of the Seven- Yea r 



Wilkinson Microwave Anisotropy Probe (WMAP7) (e.g. iJarosik et al. 



2010 



Larson et al.l |2010j). 



Here we provide an alternative estimate of the bulk using a single data set of TF measurements 
of galaxies, trimmed at faint magnitudes to ensure the linearity of the TF relation. The estimate 
is based on a method which we term ASCE for All Space Constrained Estimate. The method 
computes B(r) using ([T]) from a three dimensional t^(^) defined everywhere in a large region of 
space and constrained to match the TF data. For the analysis below, we use the SFI-I— |- surve y 



of spiral galaxies with I-band Tully-Fisher distances, (IMasters et al 



which builds on the original Spiral Field I-band Survev (iGiovanelli et al 



19991 ) and Spiral Cluster I-band Survey (IGiovanelli et al 



2006 



1997a 



Springob et al.l hoOTO 



1994 



199. 



i 



Havnes et al 



b|). We use the published SFI++ 



magnitudes and velocity widths, and derive our own peculiar velocities, rather than taking the 
published distances as given. We shall use the inverse of the Tully-Fisher (ITF) relationship. 
The main advantages of ITF methods is that samples selected by magnitude, as mos t are, will 
be rn i nimally plagued by M almquist bias effects when analyzed in the inverse direction ([Schechter 
1980l : lAaronson et al.lll982l ). We assume that the circular velocity parameter, 77 = log(Zine width), 
of a galaxy is, up to a random scatter, related to its absolute magnitude, M, by means of a linear 
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inverse Tully-Fisher (ITF) relation, i.e. 



r] = sM + rjQ. 



(2) 



The preparation of the data is done following iDavis et al.l (j201ll ). We include all field, group, 
and cluster galaxies. Galaxies in groups and clusters are treated as individual objects, though 
the redshifts for template cluster galaxies are replaced by the systematic redshift of the cluster. 
Galaxies fainter the an estimated magnitude of -20 were removed from the sample as those showed 
significant deviations from a linear TF relation. In order to get a cleaner TF sample we select only 
objects with inclination i > 45° to ease problems with inclination corrections. All this leaves us 
with a sample of 2859 galaxies with redshifts less than 100/i~^ Mpc. The effective depth of the 
sample defined as the error weighted mean redshift of galaxies is ~ 40h-^ Mpc. 

We win refer with ACD M7 to the ACDM cosmological model with the WMAP7 best fit pa- 
rameters (ILarson et al.ll201Cll ) for a flat Universe, i.e. the total mass density parameter Urn. = 0.266, 
baryonic density parameter Qi, = 0.0449, a Hubble constant h = 0.71 in units of 100 kms~^Mpc~^, 
a scalar spectral index = 0.963, and fjg = 0.8 for the rms of linear density fluctuations in 
spheres of 8h~^ Mpc. Throughout the paper, variants of ACDM7 with different ^Im and ag will be 
considered. All other parameters will be fixed at their WMAP7 values. 

The outline of the paper is as follows. Details of the ASCE method are described in ^ while 
the more standard Maximum Likelihood Estimate (MLE) outlined in ^ Tests of the methods using 
mock catalogues designed to match the SFH— |- catalogue are presented in ^ Results for the bulk 
flows from the SFI-I-+ data are given in Sj5] with the subsection ^5.11 providing a comparison with 
the ACDM models. Finally, ^discusses the results and some of their cosmological implications. 



2. The All Space Constrained Estimate (ASCE) 



Observations of distance (peculiar velocity) indicators, such as the SFI+-I- TF survey, are 
available for only a small fraction of galaxies in the local Universe (out to ~ 100/i~^Mpc). The 
absence of uniformly distributed data prevents a direct application of equation ([T]). To circumvent 
this problem, the ASCE method effectively uses ([I]) to reconstruct the bulk fiow from a 3D field 
v{x) which satisfies two conditions: a) at sufficiently large distances from the observed galaxies in 
the TF data, it has a power spectrum that is dictated by a cosmological model such as the ACDM, 
and b) it has radial peculiar velocities at the positions of the observed galaxies, which are consistent 
with the TF measurements. The approach is sini ilar to that of constrained realization from noisy 



data (Hoffman Ribak 



1991 



Zaroubi et al.lll995l ). but it is more general and easier to implement. 
Assume that the TF catalogue contains i = 1 . . . Ng galaxies with measured redshifts (in kms~^), 
czi, apparent magnitudes, rrii, and line width parameters, rji. We write the absolute magnitude of 
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a galaxy, 

Mi = Moi + P^ , (3) 

where 

Moi = mi + 51og(czi) - 15 (4) 

and 

P^ = -51og(l - Ui/cZi) (5) 

with Ui the radial peculiar velocity of the galaxy. Both czi and Ui are defined in the frame of 
the cosmic microwave background radiation (CMB). Assume that an estimate of the underlying 
cosmological velocity field, can be written as a linear combination of 

v{x) = ^a''v''{x) , (6) 

a 

where the N fields, (a = 1---N^, are gaussian random velocity fields generated using a 

cosmologically viable power spectrum. In practice these basis velocity fields will be extracted from 
a linear cosmological velocity field generated in a very large box using the power spectrum of the 
ACDM model. We then compute uf, the radial component of the at the redshift space positions 
of the observed galaxies, and define the P— basis functions, = — 51og(l — /czj) for the observed 
galaxies only. The model P is then written as 



Pt' = Y.a'^Pt. (7) 



The best fit mode coefficients a", the slope, s, and the zero point r/o, are found by minimizing the 

A'g A^a 



statistic 



\ {sMoi + sPf + m-Vif + , (8) 



i=l a=l 



where o"^ is the rms of the intrinsic scatter in 77 about the ITF relation, and TVg is the number of 
galaxies in the sample. The second term of the sum over the squares of a" is introduced in order 
to regularize the solution especially in regions of poor data coverage. In the appendix we derive 
this term from a Bayesian formulation. The solution to the equations dx^ /da°' = 0, dx^/ds = 
and dx^/drjQ = is straightforward. The coefficients a" will be used in equation ([6|) to get v{x) 
everywhere in a region of space large enough to contain the data. For each field w"(x) we compute 
its corresponding bulk flow, B"(r), according to equation ([T|) and write our ASCE bulk flow as 

BAscEW=E""B"(r). (9) 
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3. The Maximum Likelihood Estimate (MLE) 

For comparison we will present estimates of the bulk flow obtained with the standard MLE 
(jKaiserlll988i ). This method approximates the bulk flow of a sphere of radius r as the vector Bj^^e 
which renders a minimum in 

X^ = Xy^ (sMoi + 2.17s5-^ + r/o - V^) (10) 



czi<r 



with respect to the three components of B. The sum is over galaxies within r and fj is unit vector 
in the direction of galaxy i. Further, in this expression we have approximated P = — 51og(l — B • 
i-i/czi) 2.17B -r-i/czi. 



4. Tests 

In order to test the performance of ASCE and the MLE reconstructions of the bulk flow from 
the SF1++ TF data we use 2200 mock catalogs of TF measurements. In each of the catalogues, 
galaxies with the same positions as in the real SF1++ data are assigned absolute magnitudes, Mj, 
and line width parameter s, r?,;, following an artiflcial ITF relation with slope s = —0.12 and intrinsic 



scatter cx^ = 0.057 (e.g. iDavis et al.ll201ll ). The peculiar velocities of galaxies in each mock are 
taken from a linear gaussian random velocity fleld in a cubic box of 1454h"^ Mpc on the side. Each 
mock is placed randomly in this large box and the peculiar velocity of each galaxy is then obtained 
by interpolating the velocity fleld on the position of the galaxy. A gau ssian random realizaton o f 
the velocity field is generated for ACDM7 using the COSMICS package (|Ma fc BertschingeijflQoj} . 



Further, we work with a parametric form of the power spectrum taken from lEisenstein &: Hul (|l998l ) 
(eqs. 29-31 in their paper) 

For ASCE we still need to construct the basis velocity fields v°'{x). Here we use A^a = 120 
basis functions, extracted in a similar way to the mocks, but from a completely different random 
realisation of a velocity field in a very large box. Thanks to the regularization term in ([8j), the ASCE 
inferred bulk flow has very little dependence on the actual value of Ng, as long as it is large enough 
to capture the main features of the 3D flow: very similar results are obtained with N^, = 50 and 
= 120. Each of these 120 velocity flelds v°'{x) are further smoothed with a gaussian window of 
10/i~^ Mpc in width. The purpose of this small scale smoothing is to filter out low frequency modes 
which would be over-fitted by the data especially at large distances. This smoothing, however, has 
very little effect on the bulk flows reconstructed by ASCE. We emphasize that a basis function v°'{x) 
is not only defined at the galaxy positions, but also at any point in a sufficiently large volume (radius 
100h~^ Mpc) which contains all the galaxies used in the analysis. Hence, for each basis function 
we can measure its corresponding bulk fiow, B°(r), and once coefficients o" have been determined 
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Fig. 1. — Scatter plots of estimated versus true bulk flows in 400 mock catalogues. The top, 
middle and bottom panels correspond to the Supergalactic x, y and z components of the bulk 
flow. Plotted are bulk flows in a spherical region of r = 60h~^ Mpc centered at the origin. Bulk 
flows from ASCE and MLE are represented as blue dots and red crosses, respectively. The overlaid 
lines in each panel are linear regressions. 
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Fig. 2. — The same as the previous figure but where the pecuhar velocities in the mocks have been 
amphfied by a factor of 1.5. 
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e [degrees] 

Fig. 3. — The cumulative fraction of mock catalogues with estimated bulks directed within angle 
smaller than 9 from the direction of the true bulk flow. Blue solid lines and red dashed lines, 
respectively, correspond to the ASCE and MLE reconstructions. Thick and thin lines, respectively, 
refer to bulk flows of spheres of radii 40/i~^ Mpc and 100/i~^ Mpc centered on the observer. 




Fig. 4. — The differential distribution functions of the difference between estimated and true 
respective cartesian components. The notation of the lines is the same as in the previous figure. 
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from the data by minimization of ([8]) then the ASCE reconstructed bulk, B^gi^g(r), is readily given 
by (©. 

Figure [T] shows scatter plots of the estimated versus true bulk flows of a spherical region 
of 60h~^ Mpc in radius. For clarity, results from only 400 randomly selected mock catalogues 
are shown. Blue dots and red plus signs correspond to B^^e -^ascej respectively. Because 
of the anisotropic distribution of the observed galaxies, the methods may not reconstruct the 
three cartesian components equally well. Hence, x, y and z bulk flow components in Supergalactic 
coordinates are, respectively, shown in the top, middle and bottom panels. Blue and red lines in each 
panel are linear regressions of the estimated on true bulk flows. The corresponding mathematical 
expressions of the regressions are indicated in each panel. 

The regularization term in ([8]) naturally tends to underestimate the coefficients a" and sub- 
sequently the reconstructed bulk flow. However, the agreement between the ASCE reconstructed 
and the true bulk flows seen in figure ([1]) clearly demonstrates that the effect is meagre. To further 
explore the quality of the ASCE and MLE reconstructions and to ascertain that the regularization 
term does not cause a significant reduction in the amplitude of the bulk flow, we apply ASCE 
and MLE to the mock data but with true velocities amplified by a factor of 1.5. Everything else, 
including the regularization term in ([8]), remained the same. The reconstructed B^^^g and Bj^,jLg 
versus true amplified bulk flow are shown in figure ([2]). Both ASCE and MLE perform well even 
with this amplification of the bulk flow in the mocks. 

In both ASCE and MLE, the slopes of the regression lines plotted in ([1]) are close, but not 
equal to unity. The deviation from unity is significant (compared to the scatter of the points) and 
persists when the regression is done using all the 2200 mock points. This small but statistically 
significant bias depends on the radius of the sphere for which bulk flow is computed. The bias can 
easily be calibrated using the mock catalogues. Hereafter, all reconstructed bulk flows, from ASCE 
and MLE, are corrected for the systematic bias in the mean of the estimated bulk given the mean 
of the true value. In practice we write the corrected estimate of the bulk flow B^^^.^, (r) from the raw 
bulk B^^„(r) (directly reconstructed by either ASCE or MLE) as B^^j.j.(r) = CiB^^„(r) + C2 where 
Ci is the ratio of the rms values of the true to raw bulk flows and C2 is a constant term which 
accounts for the offset between the true and raw bulk flows. 

In all panels of figure ([T]), the mock B^g^.^ are tightly scattered around their corresponding 
regression lines. The scatter in B^j^^ appears to be more significant. To quantify the scatter 
between the reconstructed and true bulk flows, we plot, in figure ([3]), the cumulative fraction, 
P < {9), of mock catalogues for which the angle between estimated and true bulk flows is less than 
6. The solid blue and red dashed curves refer to ASCE and MLE, respectively, while thick and thin 
to bulk flows within 40h~^ Mpc and 100/i~^ Mpc, also respectively. The curves are computed after 
employing the correction to the systematic bias as explained above. The performance of ASCE is 
excellent. For 40/i~^ Mpc, the direction of B^g^^ is recovered within 3° for about 68% of the mocks. 
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For WOh Mpc this uncertainty increased to 7°. The ASCE method is significantly superior to 
MLE. The thin blue and thick dashed lines almost overlap, meaning that the performance of ASCE 
for lOO/i"^ Mpc is as as good as that of MLE for 40/1"^ Mpc. For r = lOO/i"^ Mpc, MLE recovers 
the direction only within 27° for 68% of the mocks. 

Figure (jH) plots the differential probability distribution function, P{5B), where 6B refers to 
the difference in all cartesian components between estimated and true bulk flows, from the 2200 
mocks. The notation of the lines is the same as in the previous figure and as displayed in the 
figure. The figure also indicates a, the rms value of 5B for the plotted cases. The low values 
of a corresponding to ASCE demonstrate its excellent ability at recovering the true bulk. The 
performance of MLE is good, but less satisfactory. 



Results 



The results of the application of the ASCE and MLE methods to recover the bulk flow from 
the real SFI++ TF catalogue are summarized in figures ([5])-(l6]). The bulk fiows are reconstructed 
for spheres centered on the Milky- Way and of radii from r = 20h~^ Mpc to 100/i~^ Mpc in steps of 
10h~^ Mpc. The smallest radius is c hosen large enough so that nonlinear effects are not expected 
to be important ()Nusser et al.lll99ll ). facilitating the comparison with cosmological models. The 



largest radius corresponds to the distance within which the data are used. Figure ([5|) shows the 
Galactic x (blue dotted), y (black solid), and z (red dot-dashed) components of B^g^^^ (top panel) 
and B^j^g (bottom) as a function of r. The magnitudes of B^gi^^ and Bj^j^^ versus radius are 
plotted, respectively, as the blue circles and plotted in figure ([6]). The la errorbars in both figures 
are based on the 2200 mocks. The component By is clearly the most significant. This is just a 
coincident and has no bearing on the statistical analysis of the results as one can always choose 
a coordinate system such that the bulk is along a given axis. The solid curve in this figure is the 
theoretical expectation of ACDM7 with but with ug = 0.85 instead of the default ug = 0.8. The 
theoretical curve is computed given the density power spectrum Pg{k, ^Im, h, Ug) by 



HIP 
27r2 



dkp,{k)W^{kr) 



(11) 



where / = with a growth index 7 ~ 0.55 for a flat Universe (jLindeii l2005l ) . and W = 
WT'B_{kr)exp{—k'^]rl?g/2) with Wth is the top-hat window function and the gaussian window takes 
care of the fact that the basis functions used in ASCE are smoothed with a gaussian win- 
dow of Rs = lO/i"^ Mpc in width. The expression (|lip is obtained assum ing the linear relation 
Hof5 = —V • v{x) between the density contrast, 5, and v{x) IPeeblesI (jl98d ). 



The MLE and ASCEreconstructed bulks are similar especially at large distances r > 40/i ^ Mpc. 
This is because, the data covers space more isotropically at larger distances. Figure ^ clearly 
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demonstrates that ASCE will not cause a significant artificial under-estimation of large bulk flows 
such as reported in lFeldman et al.l ((20101). To further, ascertain that our ASCE derived bulk flow is 
robust, the blue circles in figure ([7]) the amplitude of the ASCE bulk flow reconstructed using basis 
functions generated from a ACDM7 power spectrum but with a scalar index n = 0.75 and as = 1. 
The results are very similar to the ASCE bulk flow shown in figure ([6]) despite the significantly 
enhanced large scale power. The agreement is particularly striking at large distances. 



5.1. Comparison with cosmological models 

Figure Q indicates that the estimated bulk flows are consistent with the theoretical expecta- 
tions. But the errors are strongly correlated and a proper statistical analysis must take into account 
the covariance of the errors. Our large number (2200) of mock catalogues allows a robust deter- 
mination of the error covariance function between the bulk flow estimates at different radii. Since 
ASCE is significantly superior to MLE, we restrict the comparison with models to ASCE recon- 
struction. We use all components of B^g^^^ estimated at 8 values of r ranging from r = 30/i^^ Mpc 
to lOO/i^^ Mpc in steps of Wh^^ Mpc. The reason for not considering smaller radii is that the bulk 
is most robustly constrained independent of the assumed basis functions at r > 30/i~^ Mpc. We 
denote the set of ASCE reconstructed cartesian components at these 8 values of r by Bt and the 
corresponding underlying true quantities by Bt- We write the probability for observing the set Bq 
as 

P{Bo) = J dBtP{Bo\Bt)P{Bt) (12) 

where the probability P{Bt) for the underlying Bt is computed within the framework of a cosmo- 
logical model. Here, we adopt the ACDM model. For gaussian velocity fields, the calculation of 
P{Bi) is easily done by integrating standard analytic expressions involving the power spectrum. 
We assume that the probability P{Bo\Bt) for Bo given Bt is gaussian with error covariance matrix 
computed from the 2200 mocks. Under these assumptions, the expression (I12p yields 

where d is the number of elements in So, i.e. d = 24 = 8 x 3; for 8 values of r and 3 cartesian 
components. The d x d covariance matrix S = Sq + S^, where Sq is the covariance of the errors 
on Bo and St describes the covariance of the underlying quantities Bt- The dependence on the 
cosmological models comes through St. 
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Fig. 5. — The three Galactic cartesian components of the bulk flow as a function of radius. Top 
and bottom panels correspond to ASCE and MLE estimation, respectively. 
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Fig. 6. — The amplitude of the bulk as a function of distance for ASCE and MLE as indicated 
in the figure. The solid curve shows the rms value of the bulk flow as expected in a flat Universe 
ACDM model with flrn = 0.266, h = 0.71 and as = 0.85. 




Fig. 7. — The same as the previous figure with ASCE basis functions generated using a ACDM7 
power spectrum but with scalar index n = 0.75 and cg = 1, which has more power on large scales 
compared to our standard choice n = 0.963. 
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5.2. Consistency with ACDM7 

We begin by assessing how well ACDM7 is consistent with the data. To do that we generated 
10^ sets, Brnd, each containing d = 18 numbers selected at random from a gaussian distribution 
given by (jl3p computed with St for ACDM7. For each of those 10^ sets of B^d we compute the 
corresponding P(;Srnd) using (fT3]) and tabulate the negative of the log of the probability, nZPrnd = 
— lnP(i3rnd)- We also compute nlPo = — InP(Bo) for the observed Bo also using ACDM7. We find 
that only 26% of the 10^ values of nlP^.^^ exceed uIPq. Therefore, the ACDM7 cannot be rejected 
by the bulk flow results. 



5.3. Independent constraints on as and 7 



The ACDM expected amplitude of the bulk flow depends separately on the cosmological pa- 
ramet ers (see equation [11] and the parametric form for the power spectrum in lEisenstein fc Hu 
(jl998l )). But the most significant dependence is on fjg and and hence we restrict ourselves here 
to deriving constraints on these two parameters only. We compute nlPo for a grid of values of 
and as used in St, maintaining all other parameters at their default ACDM7 values. 

Confidence levels (CLs) on Qrn and erg are obtained by inspecting the contours of Ax^i^rm fg) = 
2(n/Po — min(n/Po)) in the (r2m,(Tg) plane. The minimum of uIPq (i.e. Ax^ = ) is at (i7m;<7g) = 
(0.236,0.88), marked by the plus sign in the figure. The ACDM7 default values (r^miCrg) = 
(0.266,0.8) are indicated by the circle. The inner and outer c ontours of Ay^ s hown in figure 
([8]) correspond to 68% and 95% CLs for two degrees or freedom (jPress et al.lll992l ). The ACDM7 
point is well within the 68% confidence level. The shape of the contours implies the correlation 
as ~ Q,^''^^. This reflects the dependence of the shape of the density power spectrum on 0,^ and 
from the factor f{^) ~ Q^'^^ (see eq. [TTI) . Only if we neglect the dependence of the shape of Pg on 
Q,m we get fjg ~ It is of interest to inspect the constraints when either of the parameters 

or cjg is fixed at certain values. Figure ([9]) shows two curves of Ax^ versus as corresponding to the 
WMAP7 ^Im = 0.266 and 0.236 giving a minimum of Ax^ in the (0m,(7g) plane as seen in figure 
([8]). Figure ([10]) plots Ax^ as a function of Qrii, for as at the WMAP7 value of 0.8 and at 0.88 
corresponding to the minimum of Ax^ in figure ([8]). In each of the curves in figures ([9]) and (|10p . 
the value of Ax^ at the minumum of the curve is set to zero. He nce, the 5x^ = 1 and 4 correspond 
to 68% (Ifj) and 95% {2a) CLs, respectively (jPress et al.l Il992l ). These curves assume a growth 
index 7 = 0.55 as is appropriate for a ACDM model. Hence figure ([9]) gives <Tg = 0.86 ± 0.11 {la) 
for ilrn = 0.266 and 7 = 0.55. However, we see from (llip that by varying 7 alone we get the scaling 
'7g(7 = 0.55) = (7g(7)^m~||'^^. W e can use this to set a constraint on 7 if we adopt fjg = 0.8 and 

= 0.266 (jLarson et al.ll20ld ). Demanding that as{'j) = 0.8, the scaling gives 7 = 0.496 ±0.096. 



Figure (llip confirms this result. The figure plots Ax^ as a function of 7 for the adopted values 



- 15 - 



of (Tg and i^rn as indicated. The left and ri ght arrows mark the v alues 7 
lower value is expected in f{R) mo dels (e.g. 



0.42 and 11/16. The 



Ga.nnouii et al.l l2009l ) and the highest corresponds to 



mouj] 
200d : 



a Dvali-Gabadadze-Porrati (DGP) (jDvali et al.l 120001 : IWeil |2008| ) flat braneworld cosmology. We 
could also substitute the scaling with 7 in the correlation ug ~ Q^'"^^ obtained from the contour 



plot to get the constraint agH.^ 



7-0.55 



(0^/0.266) 



0.28 



0.86 ± 0.11 between 7, fi^ and ag. 



6. Discussion 



The analysis presented here uses a trimmed version of the SFI++ in which galaxies fainter than 
M = —20 are removed. This ensures the linearity of the TF relation. Further, to avoid dealing with 
selecti on effects imposed on the magnitudes we use the inverse TF relation (see I Strauss &: Willick 



I995II for a thorough re view of this issue). Further, to minimize inhomogeneous Malmquist bias 



Lvnden-Bell et al.lll988l ). we do not place galaxies at their TF inferred distances, but at their 



measured redshifts which has significantly smaller observational errors. We also collapse the main 
known galaxy clusters. 

The bulk flows estimated here are remarkably featureless and do not seem to reflect the grav- 
itational effects of any of the individual main nearby clusters. Bulk flow estimates from TF-like 



relations are traditionally featureless (e.g. iDekell 1 19941 ). in contrast to velocity d ipoles estimated 



by from the distribution of galaxies in redshift surveys (e.g. iNusser &: Davislll994l ). In the analysis 
here, we have collapsed clusters and therefore, signatures of individual clusters, in our estimated 
bulk flows, could be smeared out. It is instructive to explore how much we are missing by collapsing 
clusters and wether signature of infall on clusters of nearby galaxies can actually be clearly seen 
in SFI++ or similar data. As an illustrative representative case we plot in figure [12] individual 
peculiar velocities of 54 SFI++ galaxies contained in a cylinder of 6° in radius and of 2600 kms~^ 
centered on the Virgo cluster. The individual distances, dxF to galaxies are obtained from the 
observed galaxy deviation from the straight line describing the inverse TF relation as determined 
by the 54 galaxies. The individual radial peculiar velocities are then given by Vtf = cz — HqcItf- 
Red circles in this figure show Vtf versus cz while the blue plus signs correspond to Vtf versus 
c^TF- The solid line is obtained by statistical regression of Vtf on cz, i.e. the red pints. The two 
straight dashed lines correspond to 95% confidence levels of this regression. The blue points show 
a pattern that could mistakenly be confused with actual galaxy infall onto Virgo. However, this 
pattern is entirely due to inhomogeneous Malmquist bias: Galaxies scattered to large estimated 
c^TF beyond the cluster, will also have a negative inferred Vtf- The effect of this Malmquist bias 
will be more pronounced for more distant clusters which have larger absolute errors on (iTF- The 
red points do not show a clear infall (on Virgo) of galaxies in the immediate vicinity of Virgo. 
Karachentsev &: Nasonoval (j2010l ) presents an impressive study of the observed flow of 1792 galax- 
ies near Virgo. Taken into account the inhomogeneous Malmquist bias, it is hard to detect a clear 
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Fig. 8. — Contour plot of the 68% and 95% confidence levels in the — place. The plus sign 
marks the niaximum of ttiG probs^bility distribution function at (Q^t^? (78 = (0.236,0.88), while the 
circle indicates (0.266,0.8), corresponding to the best fit WMAP7 values. 




Fig. 9. — Curves of Ax^ as a function of erg for = 0.266 (blue solid line) and = 0.236 (red 
dot-dashed) . 
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Fig. 11. — Curves of Ax^ as a function of the growth index 7 given erg = 0.8 and Q,m = 0.266. 
The left and right arrows, respectively, indicate 7 values obtained in f{R) and flat DGP models. 
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infall signature in the near vicinity of Virgo in this study as all (although the focus of their paper 
is different). 



The constraints given in figures (j8|)- (|lip show that the bulk flow alone provides useful addi- 
tional constraints on the cosmological parameters. To achieve tighter constraints one must investi- 
gate the full information in the peculiar velocity measurements. This could be done by analysis of 



power spectra and correlation functions by maximum likelihood techniques (e.g. iGorski et al 



Jaffe fc Kaiserlll995l : IZaroubi et al]ll997l : ljuszkiewicz et ahlboool : bridle et al]l200ll : lAbate fc Erdoedu 



198S 



20091 ). However, the bulk flow is particularly appealing because of its simplicity and the fact that 



it is entirely linear for sufficiently large spheres. The constraints from peculiar velocities, includ- 



probe the growth index 7 = where / is the l inear growth factor ( Peebles 



(Peebles 


1980: 


Linder 


2005) 



Adopting the WMAP7 cosmological parameters (jLarson et al.ll2010l ). we derive a local constraint 
7 = 0.495 ±0.096. This constraint is completely independent of the biasing relation between galax- 
ies and mass. Further, it is essentially a constraint at z = 0. In contrast, the lowest redshift 
constraint obtained frorn a st udy of redshift distortions in the 2dF galaxy redshift survey is at 
z ^ 0.15 (IHawkins et al.ll2003ll. Our constraint significantly improves on previous constraints on 



7 (jPossett et al.ll2010l : IWeill2008l ) derived at higher redsh ifts. This result could help us distinguish 



2008 



between alternative t heories for structure formation (e.g. lAmendola et al.ll2005l : iGuzzo et al 
Keselman et al.ll2010l ). For erg = 0.8, the constr aint on 7 disfavors PGP models at ~ 2a level, but 



it is consistent with f{R) gravity models (e.g. IStarobinskv 



2007 



Gannouii et al 



2009 



Wu et al 



2OO9I : IFu et al.l I2OIOI ) . But ug in these models should be computed self-consistently, assuming the 
same normalization at the recombination epoch. Based on WMAP7, this implies ug = 0.63 and 
0.855 for DGP and f{R) models, respectively. Adopting these erg values for these models we get 
7 = 0.315 ± 0.091 and 0.55 =b 0.098, respectively, for SGP and f (r) (Cinzia Di P orto, priv. comm). 
In the DGP model, the expected value for 7 at z = is 0.664 (IWu et al.ll2009l ). which is ruled by 
our constrain at more than the 3 a level. The f{R) model cannot be ruled out at high confidence 
level by the constraint derived here. 



Our results are in agreement with the analysis of lSandage et al.l (j2010l ) . Peculiar velo cities from 



supernovae, although very sparse, yield bulk fiows that are consistent with WMAP7 (jPai et al 



2011 



Colin et al.ll201lh . as we do. The results are in agreement with the WMAP7. The analysis of 
Bilicki et al.l (120 111 ) of the 2MASS dipole from galaxy fiuxes is also in agreement with the WMAP7 
LCPM. The direction of the bulk flow is robust and agrees with the direction of the motion of 
the local group (LG) after correctin g for the Virgq c entric infall (jSandage et al.ll2010l ). But we 
disagree strongly with the bulk flows iFeldman et al.l (j2010l ) who flnd a significant large bulk fiow 
at r = WOh~^ Mpc, using the untrimmed SFI++ survey, other individual data sets and also using 
a composite catalogue. They also use the TF estimated distances instead of the redshifts in their 
analysis of the bulk flow, leading to results which are highly susceptible to Malmquist bias. We 
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have opted to use a single uniforml y calibrated catalo g, namely the SFI++, excluding faint galaxies 



which spoil the linearity of the ITF ([Davis et al.ll201ll ) . We have also refrained from using composite 



data since minor mis-calibration errors between different catalogs could lead to large artificial flows 
when these catalogues are combined. Further, we place galaxies at their measured redshifts rather 
than estimated distances from the TF relation. This greatly suppresses inhomogeneous Malmquist 
bias which is known to lead to significant spurious signal especial at large distances. We also 
refrained from using gaussian window so that the bulk flow within a certain radius is completely 
unaffected by the increasing uncertainties at large distances. 

We have seen that the MLE and the ASCE methods give very similar results. Further, the 
ASCE bulk flow at r > 30/i~^ Mpc is almost completely independent of the cosmological model 
used in generating the basis functions. This is clearly demonstrated by the comparison of figures [6] 
andE However, in ASCE even if the results turned out to be sensitive to the assumed model used 
in generating the basis function, the validity of the model can still be confidently assessed. The 
reason is that the sensitivity would imply that the data are insufficient for constraining the bulk 
flow within the framework of the assumed model used in generating the basis function. Fortunately, 
this ambiguity is irrelevant for the SFI++ used here since the corresponding bulk flow is extremely 
insensitive to the model used in generating the basis functions. 
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A. The regularization term 

The motivation for the second term on the r.h.s in equation ([8]) can be found in a a Bayesian 
statistical formulation. Our model for the three dimensional velocity field defined at any point in 
space X is given in terms of an expansion, v^^ (x) = ^Q,a"t>"(x), over a Na basis functions v'^{x) 
each corresponding to a realization of a random gaussian field with a cosmologically motivated 
power spectrum. Given the data, the probability P(w^^|data) for v^^{x) is 

Idata) oc P{data\v^^)P{v^^) , (Al) 

where P(data|v^) is assumed gaussian with — lnP(data|u'^) = (sMoj + P/^ + r^o — 
with P^'^ related to the radial component of v^'^ by This probability function gives rise to the 
first term on the r.h.s in equation ([5D. The prior function P{v^) is the probability for the realization 
of the particular velocity field model u*^ independent of the data. For a gaussian random field 

-lnP(t;^^)= vf{x)E-\f,{x') (A2) 

X,X',J,J' 

where the indices J and J' refer to the three velocity components and J, J') is the velocity 

correlation function. Substituting the expansion = "^^=1 '^^'^^ yields 

- lnP('y^O = ^ a"a^ ^ v^{x)'~-^v^j,{x') . (A3) 

a,l3 X,X',J,J' 

Since v" are all space independent fields, the terms with a ^ /3 will be negligibly small. Since all of 
i;" are generated with the same power spectrum, the term S = J2x x' J J' ''^j{x)'^~^Vj,{x') will be 
independent of a. Therefore, — lnP(u*^) = S '^^{a'^)'^ . If the model is also to represent a random 
realization of the same power spectrum as each of the basis functions then the sum in ()A2p must 
also be equal to S. Hence P{v^^) oc P({a"}) oc exp(— X]a('^")^/2) which leads to the second term 
in (HI). 



This preprint was prepared with the A AS EflJ^X macros v5.2. 
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Fig. 12. — Individual peculiar velocities, Vtf, of galaxies in the line-of-sight to Virgo, plotted 
against the redshift (red circles) and the estimated distance ^tf (blue plus signs). The centroid of 
these galaxies is at ~ and Vtf ~ 0. 



